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We show that accounting for internal character among interacting, heterogeneous entities gener¬ 
ates rich phase transition behavior between isolation and cohesive dynamical grouping. Our analyti¬ 
cal and numerical calculations reveal different critical points arising for different character-dependent 
grouping mechanisms. These critical points move in opposite directions as the population’s diversity 
decreases. Our analytical theory helps explain why a particular class of universality is so common in 
the real world, despite fundamental differences in the underlying entities. Furthermore, it correctly 
predicts the non-monotonic temporal variation in connectivity observed recently in one such system. 


Dynamical grouping underlies myriad collective phe¬ 
nomena across the physical, biological, chemical, eco¬ 
nomic and social sciences mm- Whether the underlying 
N objects are particles, people or proteins, the issue of 
whether they evolve as isolated individuals or aggregates 
has significant consequences at the macro-level [Ml]. 
Super-radiance is driven by two-level systems coupling 
coherently via a background boson mode m-, many neu- 
rodegenerative diseases are driven by aggregation of pro¬ 
teins [Hi; large market movements are driven by traders’ 
herding |18H20j : insurgencies are driven by informal hu¬ 
man groupings mm as are gangs and online guilds 
[25j ; brain activity features collective neuronal avalanches 
[26| ; and many-body coherence phenomena are impacted 
by connectivity within exotic materials EauHi and net¬ 
works [mill]. 

It is tempting to try applying Physics models of inter¬ 
acting, identical particles to describe grouping dynam¬ 
ics in living systems. However a serious shortcoming is 
that the underlying objects (e.g. people, cells, animals 
etc.) are generally not identical, and it is this hetero¬ 
geneity that typically dictates their interactions and ul¬ 
timately their collective behavior. Even simple cells of 
a given type can have chemical, physical and conforma¬ 
tional differences that affect their interactions, while for 
humans it is usually the characteristics of other group 
members that dictate whether individuals join or leave 
a particular group [551120] • An outstanding question is 
therefore how this diversity in individual characters af¬ 
fects the dynamics of groups m la IM32] ? And how 
can this individual-level heterogeneity be reconciled with 
the emergent universality observed across many diverse 
real-world phenomena? 

This paper attempts to address these questions by 
adding a simple, continuous ‘character’ variable Xi to 
each object i, and then allowing these characters to in¬ 
fluence how objects interact with each other. We as¬ 
sign static xfs randomly from a given distribution p{x) 
though this will be generalized in future work, e.g. to 
incorporate experience. A single scalar parameter has 
already been adopted in other contexts within the so¬ 


cial science literature [H]. Here we restrict our focus to 
systems where forming links can be costly, e.g. in in¬ 
surgencies since each link increases the risk of detection 
[511 [53121], in financial trading since new links may leak 
proprietary information [181119j , in the brain since an in¬ 
crease in coordination between neurons will momentarily 
require additional resources |26j . in coherent many-body 
states in a fragmented material since each link can in¬ 
crease the chance of a decoherence event m- Likewise 
breaking a link in such systems (e.g. through a loss of 
common purpose, loss of trust, loss of coordination, or 
loss of coherence respectively) can lead to complete frag¬ 
mentation of the group (cluster) [TM551I55] - We therefore 
implement a character-driven fission-fusion mechanism 
(Fig. 1) that mimics these features, producing sparse 
networks which are visually similar to those observed em¬ 
pirically (Fig. 1(b)). Previous work [331241122] including 
in the absence of character, suggests that our main con¬ 
clusions will hold for a variety of model generalizations. 

We define the similarity between objects i and j as 
Sij = 1 — \xi — Xj\. Though we choose 0 < {xi} < 1, 
wider ranges do not affect our main conclusions. Objects 
i and j with similar characters have Sij near unity while 
those with dissimilar characters have Sij near zero. At 
each timestep t, with probability p an as-yet inexistent 
link is randomly chosen as a candidate to form. If it forms 
following the grouping rules based on Sij (see below) it 
will join together the two groups to which i and j belong 
(Fig. 1(a)). With probability 1 — p, an existing link 
is randomly chosen as a candidate to fragment. If it 
fragments following the grouping rules, the group within 
which it resides also fragments, mimicking the loss of 
common purpose, loss of trust, loss of coordination, or 
loss of coherence mentioned above (Fig. 1(a)). 

We first consider each simulation being run using one 
(and only one) of the following grouping mechanisms. 
Ml: Groups favor similar characters (e.g. kin) as in 
Fig. 1(c). At a link-forming timestep, the probabil¬ 
ity that the candidate link actually forms is Sij. At a 
link-breaking timestep, the probability that the candi¬ 
date link actually breaks is (1 — Sij). M2: Groups favor 
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much higher p. This implies that high-diversity groups 
and teams need to be encouraged by externally estab¬ 
lishing a high p (> Pc) while kinship groups will natu¬ 
rally arise for almost any p. Interestingly, the Ml and 
M2 onsets are less sharp than the intermediate M3 or 
character-free M4. This suggests that real-world popu¬ 
lations in which character dictates the grouping dynam¬ 
ics, will show far more glassy transitions indicative of 
frustrated dynamics as compared to the sharp ones in 
character-free physics models. Results shown are aver¬ 
ages over simulations with N = 10^ objects, with each 
simulation comprising 10^ timesteps and data collected 
in the steady state. The SM illustrates results for smaller 
N, and the distribution of similarities Sij for the groups 
that emerge under mechanisms Ml and M2. 

Our analytical analysis is a mean-field approach, start¬ 
ing with the coupled differential equations for Uk, the 
number of groups of size k at timestep t for k < N: 
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FIG. 1. (Color online), (a) Onr model of interacting charac¬ 
ters comprises two types of process: Link formation leading to 
joining of two groups, and link breaking leading to group frag¬ 
mentation. (b) Representative portion of PIRA insurgency 
network in Northern Ireland, adapted from Ref. [36]. Dif¬ 
ferent symbols and colors represent different character types 
(e.g. bomb-maker). It is slightly more connected than our 
model since all empirical link information is aggregated over 
a year [^. (c) Groups favoring similar characters (e.g. kin) 
illustrated by similar colors. Underneath, group size distribu¬ 
tion rik showing simulation (symbols) and analytical (lines) 
results for different p values, (d) Same as (c) but now for 
groups favoring diverse characters (e.g. team). 

diverse characters (e.g. team) as in Fig. 1(d). At a link¬ 
forming timestep, the probability that the candidate link 
actually forms is (1 — Sij). At a link-breaking timestep, 
the probability that the candidate link actually breaks is 
Sij. For comparison, we also consider intermediate (M3) 
and character-free (M4) grouping mechanisms. These are 
summarized in Table 1. 

Figure 2 shows that even for a uniform character dis¬ 
tribution p{x), rich behavior emerges. As p increases, the 
average number of links per object (A) increases from zero 
indicating groups spontaneously forming from a popu¬ 
lation of isolates. Figure 2(b) shows the corresponding 
rate of change. The position and shape of the onsets 
depend on the grouping mechanism, with the M2 onset 
(e.g. team) more abrupt than Ml (e.g. kin) but requiring 


where F is a mean-field probability of a link being formed 
between two randomly chosen objects, while Q is a mean- 
field probability that an arbitrarily chosen link will break 
and hence that group will fragment. Since our focus is 
on networks that are naturally sparse 1211 [22], we take a 
group of size k as having (/c — 1) essential links [211 [21] 
in Eq. (1), though any number 0(fc) would generate 
similar conclusions. In the steady-state, these equations 
yield two possible solutions for the number of isolated in¬ 
dividuals (see SM): m = or ni = ^■ Since 

rii < N, a transition will arise when [pF+{l—p)Q] = 2pF 
from a population comprising 100% isolates to one with 
cohesive groups, i.e. at the critical point 

Pc = QiF + Q)-^. (3) 

For p > pc^ each Uk for k > 2 changes from zero to the 
exact expression 

rife = 1^! [2fc7(fc!)(^ - k)\]-^ (dqni)''! (4) 

where 

^=pF{N-ni){N[Q{l-p)N + 2pF{N-ni)])-^ . (5) 

We can evaluate F and Q analytically to obtain pc for 
grouping mechanisms M1-M4: For a uniform character 
distribution p{x), the probability density function (PDF) 
f{y) of the similarity y = Sij is given by f{y) = 2y, with 
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FIG. 2. (Color online) (a) (A) versus p for N — 10“^ objects. 
Blue points: grouping mechanism Ml favoring similar char¬ 
acters (e.g. kinship). Orange points: grouping mechanism 
M2 favoring diverse characters (e.g. team). Purple points: 
M3 intermediate between Ml and M2. Black diamonds: M4 
character-free. Table 1 shows comparison to theoretical pc- 
(b) Rate of change, (c) pc for Ml (bottom, blue), M2 (top, 
orange), M3 and M4 (horizontal) versus inverse standard de¬ 
viation |(cr“^) of population’s character distribution p{x). 


y G [0,1]. For mechanism Ml, the probability F that 
two objects will be linked is f{y)ydy = 2/3. Similarly, 
the PDF of y associated with links is g{y) = 3t/^, hence 
the probability Q that a randomly selected link breaks 
is Jq g(y)(l — y)dy = 1/4. Following this procedure for 
M1-M4 yields the theoretical values shown in Table 1. 

Figures 1(c)-(d) and Table 1 show good agreement be¬ 
tween numerical simulation and our mean-field theory 
for {ufc} and Pc- Differences are due to neglect of higher- 
order correlations. Equation (4) further reduces (see SM) 



Pc 

FIG. 3. (Golor online) Phase diagram. Gurved phase bound¬ 
ary is our mean-field analytical result p = p^ = Q{F + Q)~^■, 
i.e. Q/F — pc(l — Pc)~^■ Diamonds show pc for uniform 
character distribution p{x). Stars show numerical results for 
gaussian p(a;) from Fig. 2(e). Ml blue, M2 orange. 


to the approximate form for p > p^'. 

k 

fc-5/2 . (6) 

Equation (6) predicts an approximately exponential cut¬ 
off at high k that depends on the grouping mechanism 
through Pc = Q{F Q)~^, together with a 5/2 power 
law exponent that does not. As data from real-world sys¬ 
tems improves, it should be possible to estimate Pc and 
p, and hence Q jF, to infer likely character-driven group¬ 
ing mechanisms in a given system. The 5/2 exponent is 
exactly that observed empirically for (i) the severity of 
attacks inflicted by insurgent groups on a civilian pop¬ 
ulation, indicating the size distribution of the insurgent 
groups [23l |24l ; (h) the distribution of stock transaction 
sizes, indicating the herd sizes of similar-minded traders 
m -, ( iii) the size distribution of neuronal avalanches, 
given avalanche initiation by a randomly chosen neuron 
(i.e. ES]); (iv) the size distribution 

of pockets of superconducting coherence in fragmented 
materials m- It is also close to the two values of 2.3 
obtained from a size study of 100 gangs in Chicago, and 
separately in Manchoukuo in 1935 m- 

Figure 2(c) shows that Pc shifts in opposite directions 
for Ml and M2 as the heterogeneity of the underlying 
population is reduced, using a gaussian character distri¬ 
bution p{x) with mean p = 0.5 and standard deviation 
a. This implies that teams require an ever higher p to 
form as a population becomes more homogeneous, with 
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Link forming 
probability 

Link breaking 
probability 

Pc (mean field theory) 

Pc (numerical) 

F 

Q 

Ml (e.g. kinship) 

S^j 

1 — Sij 

3/11 

0.10 

2/3 

1/4 

M2 (e.g. team) 

1 - s,, 

Sii 

3/5 

0.51 

1/3 

1/2 

M3 intermediate 

S^j 

Sij 

9/17 

0.49 

2/3 

3/4 

M4 character-free 

1 

1 

1/2 

0.50 

1 

1 


TABLE I. Different grouping mechanisms M1-M4. F and Q calculated analytically for uniform character distribution p{x). 
Mean-field result pc = Q{F + Q)~^. 


the population eventually comprising completely isolated 
individuals for all p. By contrast, kin groups require an 
ever lower p to form. Figure 3 shows the phase diagram. 
The numerical simulation results lie remarkably close to 
the analytic curve Q/F = pc{l —pc)~^, providing further 
support for our mean-field analysis. 

Though PIRA (Fig. 1(b)) is the best-known insur¬ 
gency network to date [36] , the data is still unfortunately 
insufficient to infer the actual grouping mechanism since 
links are aggregated over years, which is why it appears 
more dense than snapshots of our model. However we can 
test our character model against the recent state-of-the- 
art case study |36| that suggests that PIRA underwent 
a bottom-up transition over time from a rather homo¬ 
geneous organization toward team-like structures facili¬ 
tated through a process of individual contact. We start 
our model PIRA population with an MI grouping mech¬ 
anism favoring similar character links. An individual 
is introduced who uses an M2 grouping mechanism fa¬ 
voring diverse character links and hence favoring team 
formation, and who is able to spread its use to anyone 
with whom he/she instantaneously shares a cluster. They 
then become spreaders (susceptible—^infected) reflecting 
the fact that team-like structure became recognized as 
improving PIRA’s operational efficacy and hence got re¬ 
inforced over time by contacts at grass-roots level. Figure 
4 shows the resulting prediction concerning connectivity 
from our model and the actual PIRA data. The same 
non-monotonic dynamics arise in both. 

In summary, we have shown that rich phase transition 
dynamics emerge when the objects in a population pos¬ 
sess an internal character variable. Our analytical the¬ 
ory explained why a particular statistical universality is 
so ubiquitous in real-world systems, despite fundamental 
differences in the composite objects and their interac¬ 
tions. Our findings open a path toward understanding 
how different grouping mechanisms (e.g. Ml vs. M2) af¬ 
fect diseases or memes spreading in realistic (i.e. hetero¬ 
geneous) populations. In physical systems, the different 
grouping mechanisms (e.g. Ml vs. M2) can be used to 
mimic the tuning of particle-particle interactions in an 
exotic material, with p acting like an inverse tempera¬ 
ture. 

We are grateful to Chaoming Song and Stefan Wuchty 



FIG. 4. (Color online) Our model’s prediction vs. actual 
PIRA temporal variation for (top two curves) the fraction 
of isolated individuals, and (bottom two curves) the ratio be¬ 
tween the total number of network links and the total number 
of individuals. PIRA data adapted from Ref. [36] . 
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